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1 Introduction 

For a very dilute Bose-Einstein condensate which is the case in most of the 
experiments, the interatomic interaction is sufficiently weak, the mean field 
Gross-Pitaevskii theory is a logical tool to study such a system [TJ. The 
physics of the cold atom scattering is dominated by the two-body contact 
interaction which is well described by the s-wave scattering length. If we 
confine ourselves to the case of low momentum (k) and the small values of 
scattering length (a) such that ka « 1, the momentum dependence of the 
scattering amplitude of two interacting atoms can be neglected [21 [3] . How- 
ever the value of the parameter ka increases with the increase in a, even for 
small momentum (k). Recent experiments [U G2 [6] have explored the possi- 
bilities of increasing the scattering length by exploiting magnetic Feshbach 
resonance. In an experiment performed at JILA [I] it was possible to confine 
10 4 atoms of 85 Rb in a cylindrically symmetric trap at 100 nK. Around the 
Feshbach resonance at a magnetic field B ~ 155 Gauss, the scattering length 
was varied from negative to very high positive values ~ 10000 a by vary- 
ing the magnetic field. In this case the value of the parameter ka becomes 
greater than 1 at and above a = 3000 a . In this regime of ka > 1, the mo- 
mentum dependence of scattering amplitude becomes important. The first 
order theory of Beliaev [2] gives the momentum dependence of scattering 
amplitude by applying field-theoretic diagrammatic treatment to the zero- 
temperature homogeneous dilute interacting Bose gas [21 Ej. The Beliaev 
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theory was generalized to give the temperature dependence of the excitation 
spectrum of Bose gas [HI E]- In the regime ka << 1, the Beliaev's first or- 
der theory reduces to Bogoliubov spectrum [10] for contact potential. The 
work of Lee, Huang, Yang (LHY) [TTJ [T2] on the low temperature properties 
of dilute hard sphere gas considering the perturbation theory gives the first 
correction to the Bogoliubov mean field theory for dilute gases for relatively 
stronger interaction. The energy correction obtained from Beliaev's second 
order theory [2] coincides with the LHY result assuming momentum indepen- 
dent scattering length. Brueckner and Sadwa also determined the quantum 
depletion and correction to the ground state energy of a homogeneous dilute 
Bose gas [13]. To study the ground state properties of a dilute Bose gas 
with strong interactions momentum dependent scattering amplitudes along 
with the contribution from LHY correction term is important. For strongly 
interacting system (ka » 1), the higher order partial waves likely to give 
significant contribution to the scattering amplitudes. For ka > 1, only the 
even higher partial waves (e.g. d-wave ,g-wave etc.) will contribute to the 
atomic interactions in the Bose gases. It has been mentioned by Beliaev [2] 
that the d-wave contribution to energy spectrum can be ~ 10%. Here we will 
show that d-wave contribution to the column density can be double of this 
value ~ 20% at the centre of the trap for large values of scattering length. 

In this letter we present the ground state energy functional of trapped 
atomic Bose condensate including the effect of momentum dependent scat- 
tering amplitude (with full generality) and also including the LHY correction 
term. Higher partial wave interactions have been incorporated in the energy 
functional. From the ground state energy functional the time-independent 
equation has been derived to give the condensate density. We have con- 
sidered 10 4 85 Rb atoms which are cooled to a temperature of 100 nK with 
the scattering length varying from 3000 a$ to 8700 a^. The corresponding 
range of variation of ka is 1.15 to 3.33 and the values of peak gas-parameter 
x P k ~ 10~ 2 (x p k=n(0)a 3 , where n(0) is the peak density of Bose gas and a 
is the s-wave scattering length of interatomic interaction) [see Table 1]. It is 
found that for ka > 1, the contribution of the momentum dependent scatter- 
ing (considering s, d and g waves) to the column density becomes important 
and the nature of column density differs significantly from that obtained for 
k-independent scattering (contact interaction with LHY correction). The 
contribution from the g-wave interaction to the column density is negligible 
(less than 1%) for ka= 3.33. As in the vicinity of Feshbach resonance the 
two and three body loss rates play a crucial role [T5J [T5J [TBI E]> their effect 
have also been included phenomenologically in the time-independent equa- 
tions of the condensate. The patterns of the column densities as observed 
experimentally agree well with our results considering momentum dependent 
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scattering amplitude at B=157.2 Gauss (a « 3000a ). 



2 Theory 

Quantum scattering theory can be formulated in terms of partial waves. The 
wave function for relative motion is written as ip = e lkz + f(8)e lkr /r, where 
f{6) is the scattering amplitude and k is the wave vector of the scattered 
wave. The first term gives an incoming plane wave travelling along z axis 
and the second term represents a radially outgoing scattered wave [T5]. Using 
partial wave expansion the scattering amplitude for identical bosons can be 
expressed as 

oo 

f(9) = (1/k) ]T(2Z + l)exp(i5 l )stn5 l P l (cos9) (1) 
1=0 

Here, I— 0, 1, 2,... denotes the contribution of s, p, d, . . . partial waves to the 
total scattering amplitude, 5[ are the phase shifts associated with each partial 
wave, and Pi(cos6) are Legendre polynomials. For identical bosons only even 
I partial waves can contribute to the scattering amplitude. According to 
Beliaev's first order theory [2] the excitation spectrum of a dilute Bose gas 
can be given by the dispersion relation 



4 X) = V (4 + 2no/'& |) - no/(0, A* - (n f(k, A* (2) 

where e° = h 2 k 2 /(2m) is the kinetic energy and m is the atomic mass. 
Here only the forward scattering of atoms are considered. The two distinct 
scattering amplitudes (considering s-wave) are /(0, 0) = Anxa and f(k,0) = 
4irsin(ka)/k. The symmetrized amplitude f d (k/2,k/2) for higher partial 
wave (1=2) is added with f s (k/2, k/2) (for 1=0) to give the following disper- 
sion relation 



[el + 2no{(/'(£, ~) + \)} - no/ (0, 0)^P - (n f(k, 0)1)' 



(3) 

The values of symmetrized scattering amplitudes considering s and d waves 

are f s (k/2,k/2) = 47r(sin(ka)-i2sin 2 (ka/2))/k and f d (k/2, k/2) = -407r(cos5 2 sin5 2 + 

isin 2 5 2 ); where 5 2 = tan- l [{?>^cos^ - (3 - ^)sin k f }/{(3 - ^)cosf + 

3fsmf}]. 

The LHY correction term for the ground state energy of a dilute Bose 
gas [m E2] is 

r , 2tt h 2 an 128 ,na 3 ^ /9 , , 



This term can also be obtained from second order Beliaev's theory [2J. The 
ground state energy functional can be given as 

(5) 

where g x = g{2(/ s (k/2, k/2) + f d (k/2, k/2)) - /(0, 0)} and g 2 = 0); 
V tr is the trapping potential. In g x and g% the imaginary part of the scattering 
amplitudes has not been taken into account. 

By performing a functional variation with respect to ip* the Euler-Lagrange 
equation [191 EH] takes the following form 

h^V 2 + + 9M 2 + l 9 -^A\^ + ^-v^Vl 3 ]^ = M (6) 

ZiTfl 4fc 6^, O 777- 

where fi is the chemical potential which accounts for the conservation of 
number of particles. In Eq. (6) the two body and three body losses are 
introduced phenomenologically to give 

2m 4 e k 6 m 2 

(7) 

where K% and K% are the two body and three body recombination loss rate 
coefficients, respectively. 

At small momenta (ka < 1), neglecting the d-wave scattering and setting 
/ s (|, |) = f(k, 0) = /(0, 0) Eq. (6) reduces to 

+ V tr + 9l \^ + i^V^lVf ]V> = M (8) 

where g x = AnTi 2 a/m. The Eq.(8) is known as modified Gross-Pitaevskii 
equation for the atomic condensate [TP] . 

The trapping potential for the cylindrical trap with two angular frequen- 
cies u± and uj z is given as 

V tr (r) = ±m(u;lrl + uy) (9) 

The column density which is an accessible experimental quantity is de- 
fined as n c (z) = J dr±\ifj(r±, z)\ 2 . A measure of the extension of the conden- 
sate is the half width of the column density R1/2, defined as the z- value when 
n c (z = R x / 2 ) = \n c {z = 0). 



4 



Table 1: Results for the ground-state properties of 10 4 85 Rb atoms trapped in 
a cylindrically symmetric trap with ^=17.5 Hz and u z = 6.8 Hz. Chemical 
potentials are in the units of hu±; half widths are in the units of /mi. Results 
are given for k-dependent s-wave scattering [s(k)], k-dependent (s+d)-wave 
scattering [s+d(k)] and k-independent s-wave scattering [s]. 





ka 




II 

h 1 




half width 

J.K vll. VV X v_i I'll 


3000 


1.15 


s(k) 


12.47 


3.89(-3) 


21.01 






s+d(k) 


12.59 


3.84(-3) 


21.09 






s 


14.48 


3.19(-3) 


22.42 


5000 


1.92 


s(k) 


12.74 


1.49(-2) 


22.52 






s+d(k) 


13.79 


1.37(-2) 


23.09 






s 


18.77 


9.72(-3) 


25.88 


7000 


2.68 


s(k) 


11.51 


3.61(-2) 


23.72 






s+d(k) 


15.15 


2.97(-2) 


25.15 






s 


22.56 


1.96(-2) 


28.62 


8700 


3.33 


s(k) 


11.16 


5.89(-2) 


25.05 






[s+d](k) 


17.79 


4.45(-2) 


27.32 






s 


25.57 


3. 08 (-2) 


30.58 



3 Results and discussion 

Our aim here is to emphasize the significance of the momentum (k)-dependent 
scattering amplitudes (both for s-wave and higher partial wave scattering) to 
determine the ground state of the system for large scattering lengths even at 
small values of momentum. In order to do so we have considered a condensate 
of 10 4 85 Rb atoms confined in a cylindrical trap with radial (axial) frequency 
u±= 17.5 Hz (u z = 6.8 Hz) with large gas parameter values (~ 10~ 2 ) as 
achieved in the experiment of Cornish et al. [4j. We have made an attempt 
to explain some experimental results in light of the k-dependent scattering 
phenomena. In Table I we list the values of the chemical potential (/i), peak 
gas-parameter (x p k) and also the half widths of the column density distribu- 
tions (to be demonstrated in Fig. 1) along with the values of the parameter ka 
considering k-independent s-wave and also k-dependent s- and (s+d)- wave 
scattering. For the range of a considered here the parameter ka is greater 
than 1 and assumes the maximum value of 3.33 for a= 8700 a . As expected 
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the differences between k-dependent s and (s+d) results increases with the 
increase in ka and are of the order of 59% for /i, 24% for x p k and 9% for the 
half width at a = 8700ao- The differences between these results (considering 
momentum-dependent scattering) and those for k-independent s-wave scat- 
tering are also enhanced at large values of ka. Note that for the k-dependent 
s-wave scattering \i decreases with the increase in a after 5000 a . This effect 
of the negative contribution to the chemical potential is nullified due to the 
inclusion of d-wave in the k-dependent scattering amplitude. 

The column densities considering k-dependent s- and (s+d)-wave [ob- 
tained by solving Eq. (6)] and k-independent s-wave [obtained by solving Eq. 
(8)] atom-atom scattering are presented in Fig. 1(a), (b) and (c) for a— 3000 
do, 7000 a and 8700 a respectively. The dotted and solid lines correspond 
to the k-dependent s-wave and (s+d)-wave scattering results; circles are the 
k-independent s-wave results; the red dashed lines give the column densities 
with (s+d)-wave scattering considering two-body (K 2 ) and three-body (K 3 ) 
losses [obtained by solving Eq. (7)]. The values of K 2 and K 3 are taken 
from the Fig. 2 in Ref. [T5] where two- and three- body loss rates of 85 Rb 
atoms are shown as a function of magnetic field near Feshbach resonance. 
Column densities considering k-dependent s- and (s+d)-wave scattering al- 
most coincides at a = 3000a (ka = 1.15), whereas for larger values of ka, the 
(s+d)-wave scattering lowers the central column density with an expansion 
in the density distribution leading to a larger half width (as shown in Table 
1) than the column density corresponding to k-dependent s-wave scattering. 
At a = 8700 cio the deviation is 20.6% at z=0, which can be experimentally 
detected where the accuracy of measurement is higher [21]. The column 
density considering k-independent s-wave scattering is consistently lowered 
in the centre than the k-dependent results and the half width of the atomic 
cloud condensate is larger. 

In Fig. 2 we have compared our results for a = 3000 a^ with that obtained 
in the experiment [I] for B= 157.2 Gauss. In the conditions of this experi- 
ment the scattering length corresponding to B=157.2 Gauss is ~ 3000ao (as 
obtained by manual interpolation from Fig. 1 of Ref. [1]). We have solved 
Eqs. (6) and (7) by taking a = 3000a . In Fig. 2 the circles are the plots of 
the experimental data obtained by the interpolation of the condensate col- 
umn density curve taken from Fig. 3(d) in Ref. jl]. The solid line indicates 
column densities for k-dependent (s+d)- wave scattering considering two and 
three body losses [obtained by solving Eq. (7)] and the dotted line gives the 
same without considering any losses [obtained by solving Eq. (6)]. The the- 
oretical results are fitted with the experimental data at the second point 
from the z=0 axis (at z=4.03 /xra) and the deviation between experimental 
and theoretical results at z=0 is about 10%. By comparison we find that 



6 



our theoretical results are in good agreement with the experimental data up 
to z= 27 fim, but the long tail obtained in the experiment beyond 27 fim 
could not be reproduced. As the half width of the of the column density 
for k-independent s-wave interaction is larger than that for the k-dependent 
interaction [see Table 1 and Fig. 1(a)], it is obvious that the column density 
corresponding to the k-independent scattering will deviate outwardly from 
the experimental data. We have repeated the calculation with a = 3100ao 
and 2900ao (as the error in the interpolation of the curve in Fig. 1 of Ref. 
jl] is within 90 do) and found the results are almost coincident with those for 
3000 ao- Increases in the half widths obtained from our theoretical column 
densities due to the increase in a from 3000 a to 8700 a are 29.5% and 
36.4% for k-dependent (s+d)-wave and k-independent s-wave scattering [see 
Table 1]. In the experiment a= 3000 a and 8700 a correspond to B= 157.2 
Gauss and 156.4 Gauss [obtained by manual interpolation of Fig.l in Ref. 
jl]]. The increase in the half width of the column density is ~ 15% due to the 
change in B from 157.2 G to 156.4 G. Smaller increase in the half widths due 
to increase in a in the case of k-dependent results than the k-independent 
results indicates that k-independent results are closer to the experiment. In 
this study we have not considered the effect of coherent loss process and also 
the finite temperature corrections which may affect the column density for 
the large values a. 

To summarize, we have shown that the inclusion of the momentum de- 
pendent scattering amplitudes due to s-wave and higher partial waves are 
important to determine the ground state properties of the condensate when 
ka > 1 and x p k ~ 10~ 2 . Experimentally detectable significant quantitative 
differences are found between the results considering momentum independent 
and momentum dependent scattering. Theoretical column densities for 85 Rb 
atom considering (s+d)-wave scattering are found to exhibit good agreement 
with the experimental results for a = 3000 do- 

References 

[1] F. Dalfovo, S. Giorgini, L. Pitaevskii, and S. Stringari, Rev. 
Mod. Phys. 71, 463 (1999) 

[2] S. T. Beliaev, Zh. Eksp. Teor. Fiz. 34, 433 [Sov. Phys. JETP 
7, 299 (1958)] 

[3] S. Ronen, J. Phys. B: At. Mol. Opt. Phys. 42, 055301 (2009) 



7 



[4] S. L. Cornish, N. R. Claussen, J. L. Roberts, E. A. Cornell 
and C. E. Wieman, Phys. Rev. Lett. 85, 1795 (2000) 

[5] G. Roati, M. Zaccanti, C. DErrico, J. Catani, M. Modugno, 
A. Simoni, M. Inguscio, and G. Modugno, Phys. Rev. Lett. 
99, 010403 (2007) 

[6] S. E. Pollack, D. Dries, M. Junker, Y. P. Chen, T. A. Corcov- 
ilos, and R. G. Hulet, Phys. Rev. Lett. 102, 090402 (2009). 

[7] S. T. Beliaev, Zh. Eksp. Teor. Fiz. 34, 417 (1958) [Sov. Phys. 
JETP 7, 289 (1958)] 

[8] V. I. Klyatskin, Russian Phys. Journal 9, 100 (1966) 

[9] P. O. Fedichev and G. V. Shlyapnikov, Phys. Rev. Lett. 58, 
3146 (1998) 

[10] N. N. Bogoliubov, J. Phys. (Moscow) 11, 23 (1947) 

[11] T. D. Lee and C. N. Yang, Phys. Rev. 105, 1119 (1957) 

[12] T. D. Lee, K. Huang and C. N. Yang, Phys. Rev. 106, 1135 
(1957) 

[13] K. A. Brueckner and K. Sawada, Phys. Rev. 106, 1117 (1957) 

[14] H. Saito and M. Ueda, Phys. Rev. A 65, 033624 (2002) 

[15] J. L. Roberts et al., Phys. Rev. Lett. 85, 728 (2000) 

[16] P. O. Fedichev, M. W. reynolds and G. V. Shlyapnikov, Phys. 
Rev. Lett. 77, 2921 (1996) 

[17] E. Braaten and H. W. Hammer, Phys. Rev. Lett. 87, 160407 
(2001) 

[18] C. J. Pethick and H. Smith, Bose-Einstein Condensation in 
Dilute Gases, Cambridge University Press (2002) 

[19] M. Gupta and K. R. Dastidar, J. Phys. B: At. Mol. Opt. Phys. 
41, 195302 (2008) 

[20] M. Gupta and K. R. Dastidar, Phys. Rev. A 81, 063631 (2010) 

[21] D-M Stampern-Kurn et al., Phys. Rev. Lett. 81, 500 (1998) 



8 



Figure captions 



Fig. 1 Column densities of 10 4 85 Rb atoms at a = 3000 a (a), a = 7000 a 
(b) and a = 8700 ao (c) as a function of axial distance for a cylindrical 
trap with u±/2tt = 17.5 Hz and u z /2tt = 6.8 Hz. The vertical axis 

is multiplied by 2 x Oj_/N(a± = (mui±)). The dotted and solid 
lines represent the densities with k-dependent s-wave and (s+d)-wave 
scattering; circles are the k-independent s-wave results; the red dashed 
lines give the column densities considering two-body {K 2 ) and three- 
body (K 3 ) losses. 

Fig. 2 Column densities of 10 4 85 Rb atoms at a = 3000 a (corresponding 
to B=157.2 gauss in Ref. 0j) as a function of axial distance for a 
cylindrical trap with cu±/2ti = 17.5 Hz and u z /2n = 6.8 Hz. The solid 
line gives the result for k-dependent (s+d)-wave scattering considering 
two and three body losses and the dotted line gives the column density 
without considering any losses. Circles are the plots of the experimental 
condensate column density taken from Fig. 3(d) in Ref. jl]. 
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